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SUMMARY 


For  the  determination  of  coordinates  and  elevation  on  the  moon 
from  stereoscopic  photographs,  ordinary  photogrammetric  methods  are 
assumed  to  be  used.  The  highest  geometrical  quality  of  the  inter¬ 
sected  coordinates  is  to  be  expected  from  numerical  computations 
according  to  the  formulas  for  convergent  photogrammetry,  to  be  found 
in  any  of  the  more  comprehensive  textbooks  on  photogrammetry . 

The  basic  errors  of  the  procedure  emanate  from  the  image  coor¬ 
dinates  in  the  photographs  and  are  propagated  through  the  formulas 
used  according  to  the  lavs  of  error  propagation.  Regular  (systematic) 
errors  of  the  image  coordinates  are  assumed  to  be  veil  determined 
from  the  calibration  of  the  camera,  and  an  estimation  of  the 
irregular  errors  should  alvays  be  done  in  the  calibration  according 
to  the  method  of  least  squares  and  should  be  expressed  as  a  standard 
error  of  unit  veight. 

An  estimation  of  this  basic  factor  is  also  possible  from  dis¬ 
crepancies  in  the  condition  that  corresponding  reconstructed  rays 
from  the  two  photographs  shall  intersect  in  space.  An  adjustment  of 
such  discrepancies  according  to  the  method  of  least  squares  will 
give  the  standard  error  of  unit  veight  of  the  observations  behind 
the  discrepancies  in  the  intersection  condition  (in  photogrammetry, 
usually  denoted  as  y- parallaxes) .  This  standard  error  of  unit 
veight  is  closely  related  to  the  standard  error  of  unit  veight  of  the 
image  coordinates  and  can  consequently  be  used  for  an  estimation  of 
this  basic  factor  for  each  individual  pair  of  photographs. 

From  this  standard  error  of  unit  veight,  the  standard  errors 
of  all  intersected  point  coordinates  can  be  determined  according  to 
the  lavs  of  error  propagation.  In  order  to  express  these  standard 
errors  in  a  local  coordinate  system,  the  principle  of  compensation 
is  applied  to  those  intersected  points  which  are  to  be  used  for  the 
definition  of  the  local  coordinate  system.  General  expressions  for 
the  geometrical  quality  to  be  expected  in  the  planimetric  coordinates 
and  in  the  elevations  have  been  derived  and  express  this  quality  in 
terms  of  standard  errors  as  functions  of  the  basic  standard  errors 
of  unit  veight  of  the  image  coordinates. 

So  far,  no  determination  of  this  basic  factor  has  been  made  or 
reported.  The  value  used  in  the  examples,  0.01  mm,  is  estimated  with 
the  aid  of  some  preliminary  measurements  for  the  determination  of  the 
precision  of  image  coordinate  measurements  in  moon  photographs  in  a 
comparator,  i.e.,  from  repeated  and  replicated  measurements  only. 


v 


It  should  also  Toe  emphasized  that  certain  approximations, 
probably  unimportant  in  comparison  with  other  errors,  have  been  used, 
primarily  in  neglecting  the  angle  of  convergency  in  the  basic  error 
propagation  expression.  Fuxrther,  -the  coordinates  and  the  elevations 
are  referred  to  a  plane  coordinate  system.  The  transformation  into 
selenographic  coordinates  is  a  separate  operation,  not  treated  here. 

If  the  interior  orientation  of  the  photographs  is  determined 
according  to  ordinary  phot ogramme trie  methods,  the  angle  of  conver¬ 
gency  (the  libration  angle)  can  be  determined,  in  principle,  through 
the  adjustment  procedure.  The  geometrical  quality  (standard  error) 
of  this  angle  is  a  simple  fume t ion  of  the  standard  error  of  unit 
weight  of  the  image  coordinates  too.  A.  low  quality,  however,  must 
be  expected,  because  of  the  narrow  bundles  of  rays. 

The  photogrammeicric  procedure  for  the  determination  of  the 
relative  orientation  and  foar  the  estimation  of  the  geometrical  quality 
of  this  operation  can  be  directly  applied  to  moon  photographs.  If  the 
interior  orientation  of  the  moon  photographs  is  not  determined  but 
the  stereoscopic  photographs  are  taken  with  the  same  camera,  it  may 
be  possible  to  determine  the  angle  of  convergency,  provided  that 
there  are  details  in  the  plate  magazine  or  in  the  frame  of  the 
supporting  plane  whicli  can  Toe  identified  from  plate  to  plate. 

The  photogramme  trie  procedure  can  be  studied  in  available 
textbooks  together  with  the  basic  theory  of  errors  to  be  applied 
(see  for  Instance  references  [2],  [3],  and  [4]  and  the  appendix). 

Assuming  various  scales,  planimetric  and  elevation  requirements 
for  lunar  maps  and  photographs,  it  is  possible  to  formulate  a  method 
for  determining  the  dependency  of  the  geometrical  quality.  Under  the 
conditions  chosen,  it  is  shown  that;  the  lunar  instrumentation  required 
is  unrealistic.  Estimates  of  the  standard  errors  from  stereophoto  - 
grammetric  measurement  in  pti.otogra.phs  taken  with  a  telescope  of  150- 
foot  effective  focal  length  are  tabulated  for  different  values  of  the 
standard  error  of  uniic  veigtLt  (sQ)  of  image  coordinates  and  parallaxes 
on  the  scale  of  the  images.  In  addition,  estimations  of  the  relative 
geometrical  quality  f  xom  available  satellite  and  balloon  photography 
are  given. 

It  is  concluded  that  a  comprehensive  evaluation  of  the  require¬ 
ments  and  the  necessanry  geometric  quality  is  warranted  for  all  planned 
lunar  mapping  and  phoicograplay. 


GEOMETRICAL  QUALITY  OF  LUNAR  MAPPING 
BY  PHOTOGRAMMETRIC  METHODS 


I.  INTRODUCTION 

All  photogrammetric  methods  are  founded  upon  the  concept  of 
image  coordinates;  therefore,  the  geometrical  quality  of  the  final 
results  is  primarily  dependent  upon  the  geometrical  quality  of  the 
calibration  of  the  camera  and  of  the  photographs.  Cameras  for  moon 
mapping  purposes,  as  well  as  all  measurements,  must  therefore  be 
carefully  calibrated  under  operational  conditions.  In  particular, 
the  standard  error  of  unit  weight  of  the  image  coordinates  must  be 
determined  if  realistic  estimations  of  the  geometrical  quality  to  be 
expected  in  the  final  results  shall  be  made  possible.  So  far,  no 
determination  of  the  basic  standard  errors  of  unit  weight  of  the 
image  coordinates  has  been  made  or  reported. 

Stereophotogrammetric  methods  are  assumed  to  be  used  for  mapping 
the  moon  from  stereoscopic  photography  obtained  with  the  aid  of  the 
librations.  From  measurements  of  image  coordinates  and  coordinate 
differences  (parallaxes)  in  the  photographs,  ordinary  photogrammetric 
methods  can  be  used  for  the  determination  of  the  three  coordinates 
of  the  intersection  of  the  reconstructed  rays  from  the  photographs. 

The  coordinates  are  to  be  determined  primarily  according  to  the 
convergent  case  of  photography  and  are  first  obtained  in  a  system, 
the  origin  of  which  is  located  in  one  of  the  camera  stations,  usually 
the  left  one  (see  Appendix) .  The  numerical  procedure  for  the  compu¬ 
tation  of  the  coordinates  is  more  accurate  than  procedures  in  the 
ordinary  photogrammetric  plotting  instruments. 

Since  the  coordinates  are  determined  from  measured  data  via 
formula  systems,  it  is  possible  to  determine  the  geometrical  quality 
to  be  expected  through  a  study  of  the  error  propagation  from  the 
basic  measurements  according  to  usual  procedures.  The  geometrical 
quality  of  the  basic  measurements,  primarily  of  the  image  coordinates, 
must  be  determined  or  estimated  in  some  way  as  a  standard  error  of 
unit  weight.  This  fundamental  quality  can  be  determined  in  connection 
with  calibration  of ' the  camera,  but  this  calibration  must  be  made 
under  operational  conditions.  Another  possibility  is  to  use  the 
condition  that  reconstructed  pairs  of  rays  from  stereoscopic  photo¬ 
graphs  shall  intersect  in  space.  From  discrepancies  in  these  condi¬ 
tions  and  adjustments  according  to  least  squares,  using  the  elements 
of  the  relative  orientation  of  the  bundles  of  rays  and  other  possible 
"regular"  sources  of  parallaxes  as  parameters,  a  statistical  value 
of  the  residual  discrepancies  after  the  adjustment  can  be  determined. 
This  is  essentially  the  same  procedure  for  determining  the  standard 
error  of  unit  weight  of  y-parailaxes  in  aerial  photogrammetry. 
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Since  such  parallaxes  are  direct  functions  of  image  coordinates,  the  ' 
standard  error  of  the  image  coordinates  can  be  estimated  from  the 
results  of  the  adjustment  of  the  intersections  of  the  reconstructed 
pairs  of  rays.  This  method  of  determining  the  basic  accuracy  of 
image  coordinate  measurements  can  be  applied  to  arbitrary  pairs  of 
stereoscopic  survey  photographs  to  be  used  for  measuring  purposes 
and  is  very  convenient,  since  the  condition  is  furnished  by  the 
photographs  themselves  and  no  other  conditions  in  terms  of  given 
data  are  required.1  It  should  be  noted,  however,  that  there  may  be 
correlation  effects  between  the  two  bundles  of  rays  which  may  compen¬ 
sate  important  regular  (systematic)  errors  of  the  individual  rays 
and,  therefore,  the  relation  between  the  standard  errors  of  unit 
weight  of  the  image  coordinates  and  of  the  parallaxes  may  be  closer 
to  1:1  than  to  In  aerial  photogrammetry,  the  relation  1:1 

has  been  used  with  acceptable  results .  Consequently,  if  the  condition 
of  intersection  of  rays  (y-parallax  condition)  is  used  for  the  deter¬ 
mination  of  the  basic  geometrical  quality  of  the  image  coordinate 
measurements,  it  is  advisable  to  use  the  obtained  standard  error  of 
unit  weight  as  an  estimation  of  the  standard  error  of  unit  weight  of 
the  basic  image  coordinates  too.  Measurements  and  computations  for 
the  determination  of  this  basic  factor  are  comparatively  simple 
(see  references  [2]  and  [3]). 2  Stereoscopic  measuring  methods  should 
always  be  used  for  the  purpose  of  reliable  identification  in  the  two 
photographs . 

Through  differential  formulas  of  the  analytical  expressions 
used  for  the  computation  of  the  three  coordinates  of  intersected 
points  on  the  moon's  surface,  the  influence  of  the  errors  of  the 
measurements  upon  the  coordinates  can  be  studied.  If  the  errors  of 
the  measurements  are  determined  as  indicated  above,  i.e.,  as  statis¬ 
tical  expressions  (standard  errors  of  unit  weight),  the  standard 
errors  of  the  intersected  coordinates  can  be  computed  according  to 
the  special  law  of  error  propagation  (reference  [21) .  It  is  then 
assumed  that  the  errors  of  the  basic  measurements  are  normally  dis¬ 
tributed  on  a  reasonable  level,  usually  5  percent.  Tests  of  the 
statistical  distribution  of  the  residuals  in  a  sufficient  number  of 
redundant  observations  in  connection  with  the  determination  of  the 
basic  geometrical  quality  are  therefore  always  advisable.  If,  for 
some  reason,  the  basic  measurements  are  expected  to  be  of  different 
quality,  the  corresponding  (a  priori)  weights  can  be  estimated  and 
introduced  in  the  calculations.  There  is  a  clear  correlation  between 
the  resolving  power  and  the  standard  error  of  unit  weight.  The 


1The  elements  of  the  interior  orientation  must  be  known  in 
survey  photographs  and  also  the  accuracy  of  the  determination. 

2 See  also  the  Appendix.  Additional  parameters,  for  instance 
possible  "vertical"  libration,  can  be  introduced  in  the  adjustment 
procedure . 
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resolving  power  decreases  and  the  standard  error  of  unit  weight  in¬ 
creases  with  the  distance  from  the  center.  Therefore,  it  is  most 
probable  that  points  on  the  moon,  located  close  to  the  limbs,  will 
be  measured  with  less  geometrical  quality  than  those  in  the  center 
of  the  moon  and  that,  consequently,  a  weight  variation  could  be  es¬ 
timated.  Weights  can  also  be  introduced  according  to  the  intersection 
geometry.  The  size  and  directions  of  the  shadows  are  also  of  great 
importance  for  the  identification. 

The  errors  of  the  intersected  points,  determined  in  the  indi¬ 
cated  manner,  refer  to  the  origin  of  the  coordinate  system,  tem¬ 
porarily  located  in  the  left  camera  station.  For  the  determination 
of  the  geometrical  quality  of  the  mapping  procedure,  this  type  of 
errors  is  of  limited  interest  since  the  map  primarily  shall  show  the 
local  geometrical  (topographical)  relations  on  the  moon.  Therefore, 
the  local  geometrical  quality  must  be  of  primary  importance  in 
connection  with  the  mapping.  The  coordinates  and  the  elevations 
shall,  in  other  words,  be  referred  to  a  local  trigonometric  network 
or  leveling  system  on  the  moon.  Therefore,  temporary  coordinate 
systems  are  created  by  the  necessary  number  of  intersected  points, 
two  in  planimetry  and  three  in  elevation.  This  means  that  all 
intersected  points  are  transformed  numerically  into  a  coordinate 
system  defined  by  the  seven  mentioned  parameters,  which  are  repre¬ 
sented  by  suitably  located  and  well-defined  intersected  points.  The 
procedure  is  identical  with  absolute  orientation  of  a  model  in 
photogrammetry.  The  most  important  task  is  now  to  determine  the 
geometrical  quality  to  be  expected  in  the  coordinates  x  and  y  and  the 
elevations  h  in  this  local  system.  For  this  purpose,  the  principle 
of  compensation  should  be  used.  The  final  transformation  into 
selenographic  coordinates  is  not  considered  here. 

In  planning  the  lunar  mapping  project,  telescopes,  cameras, 
plotting  methods,  instruments,  operators,  and  control  density  must 
all  be  chosen  so  as  to  allow  the  coordinates  of  the  photogrammetric 
model  to  be  determined  with  a  root  mean  square  value  of  the  standard 
errors  smaller  or  equal  to  the  required  values.  General  formulas 
have  been  indicated  for  the  standard  errors  in  planimetry  and  eleva¬ 
tion.  The  standard  errors  should  be  expressed  as  functions  of  the 
standard  errors  of  the  fundamental  operations,  and  attention  should 
be  paid  to  the  balancing  effects  of  the  operations  and  of  the  control 
points.  The  number  and  location  of  the  control  points  will  influence 
the  accuracy  to  be  expected.  Therefore,  the  expressions  for  the 
accuracy  have  to  be  specialized  for  control-point  combinations. 

Wo  attention  has  been  paid  to  the  interpretation  quality 
required  of  the  photography  for  the  selection  and  interpretation  of 
the  map  details.  For  a  detailed  study  of  this  problem,  the  require¬ 
ments  on  the  map  details  must  be  specified  and  the  interpretative 
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qualities  of  the  camera  and  photographs  must  "be  investigated.  The 
resolving  power  of  the  lens -negative -diapositive -plotting  instrument- 
combination  must  be  defined  and.  tested. 


II.  APPLICATION  OF  THE  PRINCIPLE  OF  COMPENSATION  TO 
INVESTIGATIONS  OF  THE  GEOMETRICAL  QUALITY  OF 
LUNAR  MAPPING  WITH  PHOTOGRAMMETRIC  METHODS. 

1.  Determination  of  the  Accuracy  to  be  Expected  in  Elevation 
Differences  on  the  Moon.  An  assumption  is  made  that  a  photogram- 
metric  method  has  been  applied  for  the  determination  of  elevation 
differences  on  the  moon.  The  z-coordinates  (in  the  direction 
orthogonal  to  the  base)  have  been  obtained  from  the  intersection  of 
reconstructed  rays.  An  estimation  of  the  accuracy  (standard  errors) 
of  the  intersected  points  has  been  made  according  to  the  used 
numerical  formulas  and  as  a  function  of  the  standard  error  of  the 
image  coordinates.  This  standard  error  has  been  estimated  in  con¬ 
nection  with  the  calibration  of  the  camera  and/or  from  y-parallax 
determinations  and  adjustments  according  to  the  method  of  least 
squares. 

A  certain  variation  of  the  accuracy  of  intersected  points  is 
to  be  expected,  due  to  the  variations  in  the  photographic  quality 
and  identification.  The  geometrical  conditions  will  also  have 
seme  influence  upon  the  accuracy.  The  variations  of  the  accuracy 
can  be  expressed  in  terras  of  weights  to  be  assigned  to  the  z-co- 
ordinate.  In  the  following  derivation,  equal  weights  are  assumed 
for  the  sake  of  simplicity,  but  identically,  the  same  procedure 
can  be  used  for  weighted  observations  and  intersections. 

Three  elevation  points  are  selected  which  define  the  reference 
plane.  The  elevation  differences  between  the  points  can  be  assumed 
as  known,  but  with  a  degree  of  uncertainty,  since  the  points  were 
determined  from  other  measurements.  The  errors  are  to  be  of  dif¬ 
ferential  order  of  magnitude. 

The  intersected  point  group  is  transformed  to  the  elevation 
system  of  the  reference  plane.  The  errors  of  the  position  of  the 
reference  plane  will  be  entirely  compensated  in  the  three  points  by 
the  transformation  elements  which  are  one  translation,  dT|0,  and  two 
rotations,  d7]  and  dg; .  It  is  assumed  that  the  coordinates  x  and  y  of 
the  three  points  are  translated  to  a  coordinate  system,  the  origin 
of  which  is  located  at  the  center  of  gravity  of  the  three  points. 

The  differential  relation  between  elevation  errors  of  points  in  the 
plane  and  the  mentioned  parameters  is  of  the  form: 


dh  =  Xdl|  +  Yd§  +  dhQ 


h 


■where  X  and  Y  are  the  gravity  coordinates. 

Elevation  errors  dhj ,  dhg ,  and  dhg  in  the  three  points  used 
for  the  coordinate  transformation  are  related  to  the  parameters  dh^ 
dT|,  and  d§  through  the  equations  (-written  as  correction  equations): 

dT|  -  Yxdg-  dh0 

di^  =  -X^d-n  -  Yxd?-  dh0 

dhg  =  -XgdT]  -  Ysd?-  dhp 

This  system  can  be  solved  as  follows: 


dTj  [XY]  [Ydh]  -  [YY]  [Xdhl 
[XX]  [YY]  -  [XY]3 

[XY]  [Xdhl  -  [XX]  TYdh] 
d§  =  - - - - - - - 

[XX]  [YY]  -  [XY]3 


Corrections  can  now  be  computed  to  the  elevations  of  arbitrary  points 
with  the  aid  of  the  original  expression.  The  geometrical  quality  of 
the  corrections  can  be  computed  from  the  same  expression  in  applying 
the  general  law  of  error  propagation.  For  this  purpose,  the  weight 
and  correlation  number  of  the  parameters  from  the  solution  above 
have  to  be  determined.  According  to  the  definition  of  such  numbers, 
we  find: 


1 

3 


% 


m 

[XX]  [YY]  -  [XY]3 


[XX] 

[XX]  [YY]  -  [XY]3 
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M 

[XX]  [YY]  -  [XY]2 


For  points  which  are  located  in  the  vicinity  of  the  reference  plane, 
the  final  elevation  after  the  correction  can  he  expressed  hy  the 
formula: 


hfinal  =  Velimin.  +  +  xdT)  + 

Assuming  the  same  geometrical  quality  (standard  error)  in  the 
measurements  of  hpreximin  as  of  the  elevations  in  the  points  1,  2 

and  3  from  which  the  corrections  are  to  be  computed,  we  find  the 
weight  number  of  the  final  elevation  of  a  point  X. ,  Y^  in  the 
neighborhood  of  the  reference  plane: 


Si 


ihi 


[YY]  X|  +  [XX]  Y|  -  2  [XY]  X±  Y± 
[XX]  [YY]  -  [XY]2 


The  standard  error  of  the  elevation  is  then  found  from 


Bh±  *  soh  Ai.jh.j_ 

soh  is  ’the  standard  error  of  unit  weight  of  the  elevation  measure¬ 
ments  from  the  intersection  procedure.  This  error  is  primarily  a 
parallax  error  which  can  be  determined  from  the  condition  that  rays 
shall  intersect  in  space  and  then  transformed  to  an  elevation  error 
with  the  aid  of  the  distance-base  ratio.  If  the  standard  error  of 
unit  weight  of  the  parallaxes  is  denoted  s0,  the  standard  error  soh 

can  be  found  from  the  expression  sQ^  =  ^  sQ  where  d  is  the  distance 

from  the  camera  base  to  the  intersected  point  and  b  is  the  base. 

The  ratio  —  is  usually  called  the  distance-base  ratio.  For  points 
b 

located  much  above  or  below  the  reference  plane,  the  formula  system 
above  means  a  certain  approximation,  but  for  gerferal  information 
about  the  accuracy  to  be  expected  in  the  elevation  determinations, 
the  approximations  can  be  accepted.  Moreover,  the  complete  formulas 
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can  be  derived  from  the  general  rotation  or  coordinate  transformation 
formulas  without  difficulty. 

For  the  elevation  difference  between  two  points  M  and  N,  the 
weight  number  is  (with  certain  approximations  as  mentioned  above) : 


^AhAh  ^  +  ~  ^a)  +  (^m  “  ^n^2  +  “  -^n)  (^m“  Yn) 

and  the  standard  error 


SAh  "  Soh  ^QAhAh 

2.  Estimation  of  the  Standard  Error  of  Unit  Weight  and  the 
Error  Propagation.  According  to  available  information,  for  instance 
Nowicki,  reference  [1],  the  following  data  will  be  used  here: 

Average  diameter  of  the  moon  about  3, 500  km. 

Average  distance  earth-moon  about  390,000  km. 

For  the  angle  l6°  of  libration,  the  base  between  the  camera  stations 
corresponds  to  about  102,000  km.  Hence,  ^  about  3-7 

b 

Average  scale  of  available  photographs  1:22,000,000  =  1:S 
For  the  standard  error  of  unit  weight  of  parallax  measurements  s0, 
the  standard  error  of  -unit  weight  of  the  elevation  measurements  is, 
with  minor  approximation, 

soh  =  sos 


For  S  =  22  x  106  is  found 

soh  =  81  x  106  s0 


The  standard  error  of  unit  weight  of  the  parallax  measurements 
sQ,  which  is  approximately  of  the  same  magnitude  as  the  standard 
error  of  unit  weight  of  the  image  coordinates,  is  of  basic  importance, 
but  very  few,  if  any,  determinations  of  this  factor  seem  to  have 
been  made  concerning  moon  photographs.  Information  concerning  the 
geometrical  quality  of  measurements  in  star  photographs  has  been 
obtained,  and  the  figures  vary  between  2  and  5  microns.  In  moon 
photographs,  the  details  on  the  moon  are  much  less  defined  than  star 
images,  and  there  are  also  other  disturbing  factors,  such  as  different 
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shadows  in  the  two  photographs  taken  in  two  different  libration 
phases.  It  is  very  difficult  to  estimate  the  standard  error  of  unit 
weight  to  he  expected,  and  test  measurements  and  corresponding 
computations  in  accord  with  the  parallax  condition  are  highly  desir¬ 
able.  From  some  test  measurements  for  the  determination  of  the 
precision  of  image  coordinate  measurements  of  the  same  details  as 
performed  by  Nathan  Fishel,  GIMRADA,  the  standard  deviation  of  one 
measurement  was  found  to  be  of  the  order  of  magnitude  0.01  mm. 

This  figure  is  temporarily  used  as  standard  error  of  unit  weight  of 
measurements  of  image  coordinates  and  parallaxes  in  moon  photographs. 
Hence, 


soh  =  810  meters 


The  error  propagation  to  arbitrary  points  is  then  dependent 
upon  the  location  of  the  three  points  which  define  the  reference 
plane  and  the  actual  point  according  to  the  formula  given  above. 

The  lowest  standard  error  is  to  be  expected  for  Xj  =  Y.  =0.  There 
we  find 


sh  (0,0)  =  1-15soh  =  930  meters 


This  means  that  the  standard  errors  of  the  evaluations  under  the 
assumption  made  can  be  expected  to  be  of  the  order  of  magnitude 
1,000  meters,  provided  that  the  three  points  are  uniformly  distri¬ 
buted,  preferably  in  an  equilateral  triangle  of  suitable  size. 

3.  The  Coordinates  in  Planimetry.  The  differential  formulas 
of  the  coordinate  transformation  in  planimetry  are 


dx  =  dXQ  +  x  — —  ydof 

,  db 

<3y  =  dy  +  y  —  +  xdor 

°  b 


Also  in  this  case,  the  center  of  gravity  (located  halfway 
between  the  two  points  to  be  used  for  the  coordinate  transformation) 
is  used  as  origin  of  the  coordinates.  The  coordinate  errors  dXj, 
dyx ,  ,  and  dy2  of  the  two  points  are  transformed  into  the  elements 

of  the  coordinate  transformation  (absolute  orientation)  as  follows: 

dx^.  +  dx^ 

2 
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dy0 


dy^,  + 
2. 


db  =  .  +  *3^*2  +  Y1dy1  +  Y3dya 

b  Xf  +  Y®  +  Xf  +  Yf 

to  =  +  -  Vjj  -  ^ 

X®  +  Yf  +  Xf  +  Y| 


The  weight  numbers  are: 


~  2 


Qbb  “ 


b® 


Xf  +  Yf  +  X§  +  Yf 


Qto  = 


X®  +  Yf  +  X®  +  Yf 


The  influence  upon  an  arbitrary  point  of  the  errors  in  the  two 
orientation  points  can  be  found  from  an  assumed  correction  procedure. 
The  final  coordinates  of  an  arbitrary  point  X^,  Y^  can  be  found  from 
the  preliminary  coordinates  corrected  with  respect  to  the  coordinate 
transformation  in  the  two  points.  In  formulas  we  have 


-Kfinal  ~  ^prelimin.  +  +  xi  p  "  ^i^ 

Yfinal  =  Ypreiimin.  +  dy0  +  Y^  —  +  X^dh 

Assuming  all  determinations  of  the  coordinates  to  be  approximately 
of  the  same  geometrical  quality,  the  weight  numbers  of  the  final 
coordinates  are: 

3  ^  +  Yi 

‘Vi  ”  Vi  '  2  +  [XXI  ;  [YY1 
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The  standard  error  is  then  .found  from 


SX.  ~  SY±  ~  soc 

The  standard  error  of  unit  weight  is  in  this  case  assumed  to  be 
approximately  the  same  as  the  standard  error  of  unit  weight  of  the 
image  coordinates  enlarged  with  the  scale  factor. 

Assuming  the  standard  error  of  unit  weight  of  the  image 
coordinates  to  be  0.01  mm  and  the  scale  factor  22  x  10s ,  the  smal¬ 
lest  standard  error  of  any  point  is  found  for  the  point  of  gravity 
as 


sx  =  Sy  =  270  meters. 


The  irregular  errors  of  the  coordinates  can  be  expected  to  be  of  the 
order  of  magnitude  (standard  error)  of  300  meters,  provided  that  the 
distance  between  the  two  points  is  of  reasonable  length. 


III.  THE  GEOMETRICAL  QUALITY  OF  LUNAR  PHOTOGRAPHY  AND 
MAPPING  FROM  ACCURACY  REQUIREMENTS 


4.  Requirements  of  the  Lunar  Maps  According  to  Assumed  Speci- 
f  ications .  ”  . . 


Map  Scale  1: 2 50,000 
Planimetry 

Ninety  percent  of  points  that  are  to  be  depicted  on  a  map  shall  be 
located  within  800  feet  of  true  position.  No  points  shall  be  dis¬ 
placed  more  than  2,000  feet  from  true  position. 

800  feet  =  0.305  •  800  =  244  meters  (l) 

Assuming  a  normal  distribution  of  the  errors  and  a  reasonable 
level,  the  relation  between  standard  error  s  and  the  90 -percent 
condition  244  meters  is 

1.6449s  =  244 
244 

s  =  - — -  =  l48  meters  (2) 

1.6449 
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This  standard  error  is  to  be  understood  as  a  radial  standard  error. 

Assuming  the  relation  between  the  radial  standard  error  and  the 
standard  errors  of  the  x-  and  y-coordinates  to  be  J2: 1  [reference  5], 
the  standard  errors  of  these  coordinates  are 

sx  =  sy  =  105  meters  (3) 

This  standard,  error  is  to  be  understood  as  a  root  mean  square 
value  of  the  errors  of  details  on  the  map  in  relation  to  the  correct 
positions  in  a  local  coordinate  system  on  the  moon.  So  far,  this 
system  is  assumed  to  be  determined  with  the  aid  of  two  well-defined 
photogrammetrically  determined  points. 

Elevations 

From  the  condition  that  90-Percent  of  points  that  are  to  be 
indicated  on  the  moon  map  shall  be  located  within  150  feet  of  true 
vertical  position,  relative  to  an  established  lunar  datum  and  again 
assuming  a  normal  distribution  of  the  elevation  discrepancies,  the 
standard  error  of  the  elevation  is  found  to  be 


=  28  meters 


This  standard  error  is  to  be  interpreted  as  a  root  mean  square  value 
of  the  discrepancies  between  the  elevations  as  taken  from  the  map 
and  the  given  elevations  in  a  local  elevation  system  on  the  moon. 
This  latter  system  is  assumed  to  be  determined  by  some  of  the 
measured  elevations  themselves. 


Map  Scale  1: 25,000 
Planimetry 

The  90 -percent  value  is  80  feet  or  2b. b  meters.  This  cor¬ 
responds  to  a  standard  radial  error  of  l4.8  meters  and  standard 
errors  of  the  coordinates  of 

sx  =  s v  =  10 *5  meters  (5) 

Elevations 

The  90-percent  value  is  15  feet  or  b.6  meters  and  the  cor¬ 
responding  standard  error  sh  =  2.8  meters.  (6) 
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5.  Requirements  on  Geometrical  Quality  of  the  Image  Coor¬ 
dinates  .  According  to  the  derivations  made,  the  standard  errors  to 
"be  expected  in  planimetry  and  elevation  can,  with  some  generaliza¬ 
tions  and  approximations,  be  expressed  as  follows: 


SX  =  SY  ~  ^-*3s0S  (7) 

sh  =  1.2  •  3-7s0S  =  kMs0S  (8) 


s0  is  the  standard  error  of  unit  weight  of  the  image  coor¬ 
dinates  and  of  the  parallaxes.  S  is  the  scale  factor  of  the  photo¬ 
graphs,  determined  as  S  =  |  ,  where  d  is  the  distance  to  the  moon  and 
c  the  camera  constant  (focal  length  of  the  telescope) . 

Photography  for  the  Scale  1: 250,000 


Planimetry 

From  the  required  geometrical  quality  (formula  3)  and  the 
formula  for  the  quality  to  be  expected  from  the  photogrammetric 
procedure  (formula  7) ,  we  find  the  following  relations: 

1.3s0  -  =  105 
u  c 

or 

f°  =  i°L 
c  1.3d 

d  is  to  be  expressed  in  meters  and  sQ  and  c  in  the  same  units. 

Elevations 

From  formulas  k  and  8,  it  is  found  that: 

,  4.*As0  |  =  28 

£0  =  28  =  6.3063 

c  h.44d  d 

d  is  to  be  expressed  in  meters  and  s0  and  c  in  the  same  units. 


80.7692 


(9) 
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Photography  for  the  Scale  1; 25,000 


Planimetry 


From  formulas  5  and  J: 

l.3s0  §  =  10.5 

s0  _  8.07692 
T  d 

Elevations 

From  formulas  6  and  8: 

Sq  =  0.63063 

c  d 

d  is  to  he  expressed  in  meters  and  s0  and  c  in  the  same  units. 

In  all  formulas  derived,  it  is  assumed  that  convergent  photog¬ 
raphy  is  used  -when  the  lihration  is  maximum  (about  l6° ) . 

6.  Some  Applications  of  the  Formulas.  According  to  available 
information,  the  smallest  standard  error  of  unit  veight  of  image 
coordinates  and  parallaxes  to  be  expected  is  of  the  order  of  magnitude 
3  microns. 

Using  this  figure  in  formulas  (9)  and  (10)  and  the  distance  d  = 
39  ’  107  meters,  the  following  values  of  the  camera  constants  (focal 
lengths)  of  the  telescopes  -would  be  necessary: 

For  the  map  of  1:250,000: 

Planimetry 

From  (9):  c  =  14-.5  meters 

Elevations 

From  (10):  c  =  186  meters 

The  indicated  result  means  that  a  telescope  -with  a  focal  length 
of  approximately  186  meters  would  be  necessary  for  taking  the  photo¬ 
graphs  from  the  earth,  if  the  geometrical  quality  of  the  elevations 
are  to  have  a  root  mean  square  value  equal  to  28  meters.  The 


13 


telescope  dimensions  and  the  photogrammetric  mapping  of  the  moon  from 
earthbased  photography,  "with  the  mentioned  geometrical  quality,  appear 
unrealistic.  Therefore,  it  should  be  understood  that  when  the  geo¬ 
metrical  quality  is  required  to  maintain  any  standard,  the  terrestrial 
photography  intended  should  be  evaluated  comprehensively  in  order  to 
determine  -whether  the  desired  quality  can  be  obtained. 

7-  Estimation  of  the  Relative  Geometrical  Quality  in  Planimetry 
and  Elevation  on  the  Moon  from  Stereophotogrammetric  Measurements  in 
Photographs  Taken  with  a  Telescope  of  1 50 -foot  Effective  Focal  Length. 
With  some  generalizations  and  approximations,  the  standard  errors  in 
planimetry  and  elevation  on  the  moon,  after  photogrammetric  measure¬ 
ments,  can  be  expressed  as  follows: 

SX  “  Sy  1.3  5CS 

sh  =  4.44  s0S 


Here 


sQ  =  standard  error  of  unit  -weight  of  image  coordinates  and 
parallaxes,  in  the  scale  of  the  images. 

S  =  the  scale  factor  of  the  photographs. 

Maximum  libration  (about  l6°)  is  assumed,  and  further,  it  is 
assumed  that  the  errors  of  the  image  coordinates  and  parallaxes  are 
of  irregular  nature  and  are  normally  distributed  at  a  reasonable 
level. 


The  scale  factor  S  can  be  expressed  as  S  =  -  where  d  is  the 
distance  to  the  moon  and  c  is  the  effective  focal clength  (camera 
constant)  of  the  telescope. 

d  is  of  the  order  of  magnitude  39  *  107  meters,  and  c  is 
assumed  to  be  150  feet,  or  about  50  meters . 

Hence,  the  standard  errors  in  planimetry  and  elevation  can  be 
written: 


1.3sn  •  39  ’  107  „ 

sx  =  Sy  =  - 2 - - -  =  1.01  •  107  S0 

50 


s .  = 


4.44s.  •  39  •  10 


=  3*46  •  l07s. 


50 
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The  standard  errors  are  tabulated  for  different  values  of  sQ. 
Standard  errors  in  planimetry  and  elevation  for  different  values  of 
s0  are: 


so 

(mm) 

sx 

SY 

(meter) 

Sh 

(meter) 

0.001 

10 

35 

0.002 

20 

70 

0.003 

30 

io4* 

0.004 

4o 

138 

0.005 

50 

173 

0.006 

6o 

208 

0.007 

70 

242 

0.008 

80 

277 

0.009 

90 

311 

0.010 

100 

346 

0.015 

150 

519 

0.020 

200 

692 

0.030 

300 

1038 

*This  value  of  sfl  is  the  lowest  to  be 
expected . 


From  calibration  of  the  telescope  and  of  the  photographs  under 
operational  conditions,  the  sQ  value  for  the  instrument  is  to  be 
determined.  From  parallax  measurements  in  stereoscopic  photographs, 
which  are  to  be  used  for  the  mapping,  the  actual  value  of  s0  can  be 
estimated  according  to  the  method  of  least  squares.  Confidence 
intervals  can  be  estimated  from  the  number  of  redundant  observations 
(degrees  of  freedom)  and  the  confidence  level,  usually  5  percent. 

8.  Estimations  of  the  Relative  Geometrical  Quality  from 
Available  Satellite  and  Balloon  Photography.  The  highest  geometrical 
quality  to  be  expected  in  lunar  mapping  from  available  photographs 
can  be  estimated  from  the  formulas  derived. 

Assuming  the  scale  1:22,000,000  in  the  photographs  and  s0  = 
0.003  mm.,  the  formulas  (7)  and  (8)  give  the  order  of  magnitude  of 
the  standard  errors  to  be  expected,  as  follows: 

Sx  =  Sy  =  86  meters  =  282  feet 
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=  293  meters  =  960  feet 


For  photography  from  satellites  or  balloons,  similar  principles 
can  be  used,  as  shown  above.  In  order  to  fulfill  minimum  elevation 
requirements  for  lunar  mapping,  stereoscopic  photography  from  satel¬ 
lites  or  lunar  orbitors  appears  to  be  the  best  method.  With  reference 
to  the  plotting  of  the  photographs,  the  actual  geometrical  quality 
to  be  expected  can  be  estimated  [Section  I]  from  redundant  parallax 
observations  and  computations.  If  "vertical"  stereoscopic  photography 
is  to  be  used,  the  basic  formulas  must  be  modified. 

Finally,  it  should  be  emphasized  again,  if  suitable  lunar  maps 
are  to  be  made,  then  calibration  of  the  cameras  and  tests  of  the 
geometrical  quality  of  image  coordinates  of  the  lunar  photographs 
must  be  made  under  operational  conditions  with  clearly  defined 
requirements  for  the  geometrical  quality  that  must  be  fulfilled. 
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APPENDIX 


ELEMENTARY  DERIVATION  OF  FUNDAMENTAL  FORMULA  SYSTEMS  FOR  ANALYTICAL 

PHOTOGRAMMETRY 


1.  Independent  pairs  of  photographs. 


It  is  temporarily  assumed  that  the  normal  case  of  photogrammetry 
is  strictly  present,  i.e. ,  the  camera  axes  are  parallel  and  orthogonal 
to  the  base,  Fig.  1.  The  image  coordinates  are,  for  reasons  which 
will  be  explained  later  on,  written  with  parentheses  (x1),  (y1)  in  the 
left  and  (x"),  (y")  in  the  right  photograph.  The  origins  of  the 
image  coordinate  systems  are  located  in  the  principal  points  (H')  and 
(H") ,  respectively,  and  the  principal  distances  (camera  constants  or 
calibrated  focal  lengths)  are  denoted  (a,)  and  (ca),  respectively. 

The  central  projections  are  temporarily  assumed  to  be  mathematically 
correct,  and  all  other  data  are  assumed  to  be  errorless. 

Because  the  principal  distances  (cx )  and  (cg )  are  assumed  to  be 
different,  the  image  coordinates  of  one  of  the  photographs,  for 
instance  those  of  the  right  photograph,  are  transformed  to  a  photo¬ 
graph  with  the  principal  distance  (c^ )  as  in  the  left  photograph  of 
Fig.  1.  The  transformation  is  a  simple  similarity  transformation 
and  is  made  through  multiplying  the  image  coordinates  (x"),  (y")  by 
the  factor 


hi 

(°a) 

The  transformed  coordinates  are  then  transferred  to  the  left 
photograph  as  shown  in  Fig.  1.  Then,  immediately  obtained  from 
similarity  is 


_ 

(x')  -  J£Lis.> 


(ci)  (c5) 

(*’)  (c2)  -  (x”)  (Cl) 


(1) 
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Fig.  1.  The  normal  case.  Independent  pair  of  photographs.  Deriva¬ 
tion  of  the  relation  between  transformed  image  coordinates  (x')> 
(y')>  (x")>  (y")  and-  model  coordinates  x,  j,  z.  The  image  coor¬ 
dinates  of  the  right  photograph  are  first  transformed  to  the  prin¬ 
cipal  distance  (c^ )  of  the  left  photograph.  Then  the  fundamental 
principles  of  the  normal  case  of  stereophotogrammetry  are  applied. 
Note  that  the  principal  distances  (Cj)  and  (c2)  as  well  as  the 
positions  of  the  principal  points  (H')  and  (H")  are  functions  of  the 
image  coordinates  and  the  elements  of  the  interior  and  relative 
orientation  of  the  rotated  photographs. 


X 


(2) 


(x')z  _ 

(c3)  (x') 

(cx) 

(x*)  (c2)-(x")  (cx) 

(y')z 

(s)  (y) 

(Cl) 

'  (x*)  (c2)-(x")  (cj 

(y")z 

K)  (y") 

(s) 

(x')  (c2)-(x")  (CX) 

(3) 

W 


These  coordinates  are  evidently  determined  in  a  local  coordinate 
system,  the  origin  of  'which  is  located  in  the  left  exposure  station 
0X .  Through  coordinate  transformation  they  can  be  transformed  and 
expressed  in  another  coordinate  system  (absolute  orientation)  with 
the  aid  of  a  sufficient  number  of  control  points  or  other  control 
data,  at  least  seven  parameters.  The  condition  that  corresponding 
rays  must  intersect  can  evidently  be  written 

yx  =  y3  (5) 

or  according  to  (3)  and  (4) 


(y')  -  jyo  =  o 
<°i)  (s> 

(y’)  (c2)  =  (y”)  (cx) 


(6) 


(7) 


If  two  rays  do  not  intersect,  a  y-parallax  (vertical  parallax 
according  to  terrestrial  photogrammetry)  will  appear,  defined  as 

py  =  yi  -  y3  (8) 
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After  substitution  of  the  expressions  (3)  and  (4)  is  found 


„  -  r  (y’_)  _  (y")  1 (y*)  (*3)  -  (y")  (ci)  >, 

y  w  (s) i  ' srsr 


(9) 


If  the  parallax  is  assumed  to  be  caused  by  small  errors,  the 
following  expression  is  generally  used 


py  =  dyi  -  dya 

1.1  Determination  of  the  relative  orientation. 


If  the  elements  of  the  relative  orientation  of  the  photographs 
are  not  sufficiently  known  or  are  completely  unknown,  the  expression 
(9)  can  be  used  for  a  determination  of  the  elements.  There  are 
five  elements  of  the  relative  orientation,  and  relating  to  these 
elements  are  assumed  conditions  requiring  five  angles (independent 
pairs  of  photographs) .  Both  the  exposure  stations  are,  in  other 
words,  assumed  to  be  fixed. 

The  five  angles  are  usually  chosen:  two  angles  (q\  and  k1) 
of  the  left  and  three  angles  (y^ ,  c^,  ^ )  of  the  right  photograph. 

If  a  photograph  in  the  moment  of  exposure  was  rotated  through  three 
angles  around  perpendicular  axes  through  the  perspective  center 
(Fig.  2)  and  the  image  coordinates  Xy,  yr  and  the  principal  distance 
c  are  given  for  a  point  P  in  the  rotated  photograph,  the  correspond¬ 
ing  coordinates  (xr),  (yr) ,  (c)  of  point  P  can  be  expressed  in  the 
original  coordinate  system  (before  the  rotation)  with  the  aid  of  a 
well-known  coordinate  transformation  through  rotation  of  axes.  The 
actual  transformation  fozmulas  depend  upon  the  order  of  the  rota¬ 
tions  and  the  sign  of  the  angles .  Here  those  arrangements  and  signs 
will  be  used  which  are  to  be  found  in  the  Wild  Autographs  A7  through 
A9  (uj  primary,  9  secondary,  h  tertiary),  but  completely  arbitrary 
arrangements  and  signs  can  be  assumed.  The  transformation  formulas 
can  be  derived  with  different  methods,  see  for  instance  reference 
[1]  where  a  very  elementary  but  completely  general  procedure  has 
been  applied. 

For  the  conditions  as  shown  in  Figs.  2  and  3,  the  transforma¬ 
tion  formulas  from  the  system  xr,  yr,  c  into  the  system  (xr) ,  (c) 
are : 
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Fig.  3.  Transformation  of  the  image  coordinates  x',  y' ,  c  of  the 
left  photograph  to  the  coordinates  (x')>  (y')>  (c)  in  the  x-,  y- 
z-system. 


(Xy)  =  Xj.  coscpcosK  -  yrcoscpsinK  -  c  sincp  (lOa) 

(y  )  =  xr  (sinmsinTcosx  +  costusinn)  + 

+  yr  (cosujcosk  -  sinuasincpsinx)  +  c  sinoocoscp  (10b) 


(c)  =  xr  (-sinuisinn  +  costusincpcosx.)  - 

-  yr  (cosuosincpsinH  +  sinoocosn)  +  c  cosuocoscp  (10c) 

It  should  be  noted  that  the  coordinate  (c)  will  vary  from  point 
to  point  and  that,  consequently,  if  the  transformation  is  regarded 
as  a  transformation  of  image  coordinates  from  a  rotated  photograph 
to  a  not  rotated  photograph,  this  latter  one  will  obtain  a  specific 
principal  distance  (c)  for  each  individual  point,  Fig.  3- 

The  expressions  (10a)  through  (10c)  are  next  applied  to  the 
two  photographs,  the  angles  of  rotations  of  which  are  cpx  ,  ^  and 
U)s,  cpa ,  Kg,  respectively,  cox  can  be  chosen  arbitrarily  because  so 
far  only  the  relative  orientation  is  considered.  It  can  therefore 
be  chosen  zero,  and  u)3  is  consequently  temporarily  to  be  regarded 
as  a  difference  only. 

From  measured  image  coordinates  x',  y*  and  x'1 ,  y"  and  the 
principal  distance  c  (equal  for  the  two  photographs)  are  obtained 


(x>) 

=  x’cosTjCoskj  - 

y  'coscf^  sinKx  -  c 

sin9x 

(lla) 

(y) 

=  x’sinKj  +  y'coSHj 

(Hb) 

(cj 

=  x'sin^cosKj  - 

y'sincpjSin^  +  c 

coscpj 

(lie) 

(x") 

=  x"coscp2cos»tj  - 

y"cos9gSinKg  -  c 

sin9a 

(12a) 

(y:) 

=  x"  (sinm2  sintfe  cos^  +  coso^sinKg) 

+ 

+  y'^cosuigCosKg  • 

-  sinocfeSincf2sinKg) 

+  c  sinojgcoscffe 

(12b) 

(cs)  =  x'^-sinu^sinKg  +  cosujjSincDgCoSMg ) 
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-  y"(cosw2sin^2  sinxg  +  sinu>2cos?<g )  +  c  cosuu2  coscps  (12c) 


If  the  angles  of  the  relative  orientation  are  assumed  to  be 
known  with  sufficient  quality,  the  expressions  (ll)  through  (l2)  can 
directly  be  computed  for  each  point  and  then  substituted  into  the 
expressions  (l)  through  (4).  The  conditions  (5)  through  (7)  can  be 
applied  to  every  pair  of  rays  and  serve  as  checks  that  angles,  image 
coordinates,  frequently  other  data,  and  the  calculations  are  correct. 

Because  there  are  always  inevitable  sources  of  errors,  certain 
discrepancies  must  always  be  expected  in  the  condition  of  intersec¬ 
tion  of  pairs  of  rays.  Such  discrepancies  can  then  be  used  for  the 
determination  of  the  sources  of  errors.  Usually,  and  in  particular 
in  aerial  photogremmetry  but  sometimes  in  terrestrial  photogrammetry 
too,  the  angles  used  in  the  first  calculation  are  approximate  and 
are  consequently  regarded  to  be  affected  with  errors,  which  must 
be  determined.  First  the  discrepancy  p  is  computed  in  at  least  five 
suitably  located  points  with  the  expression  (9) .  These  discrepancies 
are  then  used  for  the  determination  of  corrections  to  the  approximate 
angles.  Since  the  functions  to  be  used  -  the  expression  (9)  after 
substitution  of  the  expressions  (lla)  through  (12c)  -  are  trigonome¬ 
trical  (transcendent),  a  direct  solution  is  generally  not  possible. 

An  iterative  procedure  must  generally  be  used.  First,  differential 
formulas  of  the  complete  expression  (9)  are  developed  around  the 
approximate  values  of  the  angles.  Then  the  necessary  five  linear 
equations  can  be  set  up  and  solved,  provided  that  the  general  condi¬ 
tions  for  the  solution  of  such  an  equation  system  are  fulfilled. 

If  the  found  corrections  to  the  preliminary  values  of  the  angles  are 
too  large,  the  procedure  is  repeated  after  a  differentiation  around 
the  corrected  values  of  the  angles.  The  procedure  is  iterated  until 
the  corrections  are  smaller  than  certain  tolerances,  which  can  be 
determined  with  the  aid  of  statistical  principles  and  with  respect 
to  the  basic  accuracy  of  the  procedure.  This  accuracy  is  to  be 
determined  from  redundant  observations  and  with  the  aid  of  the  method 
of  least  squares. 

As  an  example,  the  derivation  of  linear  differential  formulas 
of  the  expression  (9)  will  next  be  shown  for  the  important  normal 
case,  i.e.,  under  the  assumption  that  the  approximate  values  of  the 
angles  are  zero.  Under  such  assumptions,  cosines  and  sines  in  the 
expressions  (ll)  through  (12)  are  to  be  substituted  by  1  and  the 
differentials  of  the  angles,  respectively.  Hence 


(x1)  =  x'  -  y'dHx  -  cd9x 


(13a) 
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(y')  =  y'  +  x’d^  (13b) 

(CX)  =  C  +  x '  dcpx  (13c) 

(x")  =  x"-  y"dMg-  cd cp2  (l4a) 

(y")  =  y"+  x"dng  +  cdu^  (l4b) 

(c2)  =  c+  x"dcp3  -  y"du)s  ( l4c ) 


All  squares  and  products  of  differentials  are  neglected. 
After  substitution  of  (13)  and  (ik)  into  (9)  is  obtained 


y!-y2+Py  = 


x'y" 

y,-y"+x'dx1 - - —  dcpx  -  x"dHg  + 


x"y' 


y'v"  > 

dcfe-  (i+  2—r-  )  CdOL 


r-v-tt 


-y'd^-  (1+  )  cdcp1+y"dH3  +  (l+  3L^_)cdc^ 


(x'-x")  (1  + 


x’y" 


c3  /  — U  ~"3  c3  '-—rg  c 


d(J-b 


(15) 


Since  the  fraction  in  the  denominator  under  normal  circum¬ 
stances  is  smaller  than  1,  the  entire  expression  can  be  developed 
into  a  geometrical  series,  according  to  the  well-known  expression 


S  = 


a 


1-9 


=  a  +  aq  +  aqs 

(9  <  l) 


+ 


Terms  containing  differentials  of  second  or  higher  order  are 
neglected.  After  simple  developments  and  -with  minor  approximations, 
the  following  well-known  y-parallax  formula  is  found  from  (15) 


x'-x" 


=  p*  =x,d*1-  Xc^  dcp1-x"dK3+  dcp2-(l+  y  y-) cd<A 


(16) 


It  should  be  noted  that  x'-x"  approximately  is  the  base  b'  in 
the  photographs  for  limited  elevation  (depth)  differences  in  the 

b* 

object  and  that,  consequently,  the  fraction  ^  approximately 
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C  JL 

represents  the  scale  —  =  —  of  the  photographs.  Further  should  he 

Z  D 

noted  that  p  in  (15)  represents  that  part  of  the  total  y-parallax, 
which  is  caused  by  the  approximation  of  the  angles  only  and  that  it 

X  *  —  xn 

should  be  defined  as  py  =  y:  -y2  .  The  expression  Py  — — —  in  (l6)  is 

the  parallax  according  to  (9),  transformed  to  the  image  scale  "where 
it  is  denoted  p'. 

For  the  approximate  normal  case  and  limited  elevation  (depth) 
differences,  certain  simplifications  can  be  introduced  in  (16) . 
Assuming  y'=y"  and  x"=  x'-  b'  is  found 


y'-y"  =  p'  =  x'd^-  X(~~  dc^-  (x'~  b')dM2+  —  —  dcp3 -  (l+^ 

y  ^  ^ 


-)  cdtOg 

(17) 


These  differential  formulas  are  very  well  known  in  photogram- 
metry  and  must  be  regarded  as  the  most  important  mathematical  expres¬ 
sions  for  the  important  normal  case  in  practice.  The  final  values 
of  the  elements  of  the  relative  orientation  are  to  be  determined 
indirectly  from  the  basic  observations  of  image  coordinates,  with 
the  aid  of  the  elements  of  the  interior  orientation  and  using  the 
condition  that  corresponding  rays  shall  intersect.  This  procedure 
is  identical  with  the  practical  procedure  of  the  relative  orientation 
in  a  stereoscopic  restitution  instrument.  The  accuracy  of  the  pro¬ 
cedure  can  uniquely  be  determined  from  redundant  observations  accord¬ 
ing  to  the  method  of  least  squares,  and  tolerances  for  the  residuals 
of  the  elements  of  the  relative  orientation  and  other  data  can  be 
determined  accordingly.  The  elements  of  the  interior  orientation, 
including  possible  regular  errors  of  the  central  projection  (image 
coordinates)  and  the  irregular  errors  of  the  image  coordinates  are 
evidently  of  basic  importance  for  the  accuracy  of  the  determination 
of  the  elements  of  the  relative  orientation  and  all  functions  of 
these  elements.  The  influence  of  errors  of  the  elements  of  the 
interior  orientation  upon  the  elements  of  the  relative  orientation 
can,  for  instance,  be  found  from  differentiations  of  the  complete 
analytical  expressions  (ll) ,  (12),  (9)  and  then  application  of  the 
laws  of  error  propagation.  It  should  be  noted  that  the  image  co¬ 
ordinates  x',y*  and  x",y"  in  the  expressions  (lla)  through  (l2c)  in 
principle  are  differences  between  measured  image  coordinates  (sub¬ 
script  m)  and  the  coordinates  of  the  principal  points  (subscript  o), 
or  in  formulas 

x'  =  *£  -  x"  =  x^  -  x£ 
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y  =  k  -  ^ 


V  =  V  '  -  v" 
J  Jm  Jo 


Consequently,  the  propagation  of  errors  in  the  coordinates  of 
the  principal  points  can  he  studied  from  these  expressions  and  dif¬ 
ferentiations.  It  must  be  emphasized  that  the  errors  of  the  prin¬ 
cipal  distance  and  of  the  principal  points  as  well  as  of  other  ele¬ 
ments  of  the  interior  orientation  (distortions,  etc.)  must  be  deter¬ 
mined  in  connection  with  the  calibration  of  the  camera  and  of  the 
photographs ,  which  must  be  performed  under  real  operational  condi¬ 
tions,  in  addition  to  laboratory  tests.  Tolerances  must  also  be 
specified. 

The  relative  orientation  and  the  condition  that  corresponding 
pairs  of  rays  shall  intersect,  (5)  through  (9),  is  evidently  a  very 
valuable  check  of  certain  regular  errors  of  the  elements  of  the 
interior  orientation  and  also  makes  possible  a  determination  of 
certain  regular  non-pro jective  errors,  which  are  not  correlated 
between  corresponding  points  in  the  two  photographs.  The  residual 
irregular  discrepancies  of  the  intersection  condition  can  finally  be 
estimated  as  a  standard  error  of  unit  weight  and  can  frequently  be 
used  for  investigations  of  the  error  propagation  in  functions  of  the 
elements  of  the  relative  orientation,  i.e. ,  the  results  of  the  entire 
photogrammetric  procedure. 


1.2  Determination  of  model  coordinates. 


After  the  determination  of  the  elements  of  the  relative 
orientation  according  to  the  iterative  procedure  as  described  above, 
the  coordinates  z,  x,  and  y  can  be  computed  on  a  temporary  scale 
according  to  (l)  through  (4) .  These  coordinates  are  termed  model 
coordinates .  Because  the  two  y-coordinates  yx  and  y3  generally  are 
somewhat  different  due  to  residual  y -parallaxes,  the  average  of  the 
two  values  is  used  as  the  final  y-coordinate,  i.e. 


y  = 


y  +  y3 
2 


(18) 


It  would  be  more  correct  to  determine  weights  of  y1  and  y3 
according  to  (3)  and  (4)  and  then  to  determine  the  weighted  average. 

For  many  purposes,  the  differential  formulas  of  (l)  through 
(4)  are  of  great  interest.  The  formulas  can  be  derived  in  the  same 
manner  as  was  shown  above  concerning  the  y-parallaxes . 
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Hence 


dz  -  z  d*'+  (1  +  W^)  ** -  z  - 

-  (1+  £i)  £  -  zdtp  +  _££_  zd^  (19) 

C  (x'-x"  2  v"V  ^ 


.12 


(x' -x")c 


It  should  be  noted  that  z  is  positive  in  the  direction  of  the 
photography.  The  elevation  error  dh  in  a  model  of  approximate  ver¬ 
tical  photographs  must  therefore  have  reversed  sign  in  comparison 
■with  (19)  if  the  elevations  are  positive  upwards. 

Further  are  found 


dx  =  z 


x  y  vi2  .y.«t  x'v” 

(*'-*")  <=  dKl+  (1  +  ^  (*'-*")  Zdtfl'  Z  (x'-x")  c  d,<3  ' 

-  u  -  ?)  ^  (20) 


c2  (x'-x") 


(x'-x")c2 


r  I  2 


dy  =  (x'  +  y'--0  |  d^-h  (1  +  ~~)—L l5_  dtp,  -  /  y  ?  d Hr, 

1  x'-x"  c  ^  c2  x'-x"  Yl  (x'-x  )c  ^ 


y'y"z 


.  (1  +  2e2)  r±.  d9a+ 

-*2  X*-X”  3 


(x*-x")c2 


(21) 


dy3,  dK,+  (1  +  5l)  -£*  -  dcp,-  (-JCL  _  X")  «  dH, 

3  (x'-x")c  ^  c2  X'-x"  T1  x'-x"  c  * 


(1  +  ^f)  4^7  d^+  ^  +  (22) 

c2  x'-x  (x.'-x")  c2 


If  the  final  y-coordinate  is  determined  according  to  (18) ,  the 
corresponding  differential  formula  is 


dy  =  .t_dy£ 
2 


(23) 
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2.  Dependent  pairs  of  photographs  (extension) . 


In  this  case,  the  elements  of  the  relative  orientation  of  one 
(usually  the  right)  photograph  only  are  used.  Consequently,  the 
two  translations  by  and  bz  must  also  be  used  together  with  the  angles 

<Hj  ,  %,  and  *8  • 

According  to  Fig.  4,  the  base  b  is  divided  into  three  components 
bx,  by,  and  bz.  The  directions  and  the  arrangement  of  the  transla¬ 
tions  are  chosen  according  to  the  Wild  Autographs. 

Directly  from  Fig.  4,  the  model  coordinates  can  be  determined 
from  similarity.  The  base  component  bx  must  be  corrected  with  the 

term  -  (x")^  ,  which  is  also  derived  from  similarity. 

(Gs) 

The  following  expressions  are  found 


(x'QbZa 

(eg)  _  ( Cg ) bx  -  (x")bzg 

(x")cx  1  (Cg)x'  -  (x")  cx 

(C3) 


_  x'z  _  (ca)bx  -  (x'QbZg  x, 
Cj  <c3)x'-  (x'Ocj 


y'z  _  (c5)bx  -  (x")bzg 

ci  (cs)x'  "  (x")  ci 


bx  - 
z  =  - 

X'  - 


(30) 


(3D 


(32) 


y2= 


(y")(z  -  bZg  ) 

”  (Cg) 


by8 


cxbx  -  x'bZg 
(c3)x'-(x")c1 


(y") 


(33) 


The  expressions  for  (x"),  (y"),  and  (cg)  are  given  in  (l2a) 
through  (12c) .  If  the  same  camera  was  used  for  taking  both  photo¬ 
graphs,  cx  =  c. 
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Fig.  Ij-.  The  normal  case.  Extension.  The  influence  of  the  "base 
component  bZg . 
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Further  is  found 


(y")(z-hz^) 

M 


+  by2 


(34) 


After  substitution  of  the  expressions  (32)  and  (33)  is  found 


P  = 


bx[(c2)y'-(y")c1]  4-bZj,  {(y")x«-(x")y^  +by2{(c2)x'-(x”)c1  ] 

(cs)x>  -  (x")^ 


(35) 


Differential  formulas  can  be  derived  as  shown  above,  see  (15) 
through  (17) .  For  the  approximate  normal  case  and  limited  elevation 
(depth)  differences  in  the  object  is  found 


Py  =  Py  =  -x"d*3  +  2L2L  dcp2  -  (1  +  2L_  )cdco3 


+  —  dbz2+  |  d^ya  (36) 


Further  are 


-rtt  J.V*  Y**  rr  Y** 

dz  =  -  ^  z^+(x"-  >  5  ^  dbz*+ 


x"y" 


(x'-x")  ^(37) 


Note  that  dz  =  -dh  for  approximately  vertical  aerial  photo¬ 
graphs  . 


dx 


x'y" 

(x'-x”) 


z  ,  (x,,2-i-  c2)b 

d i x'dcp3- 


(x'-x")2  c 


Jg—  dbZg  +  x'x"y1-z  „dur 
(x'-x")c  ^  (x'-x")^  3 

(38) 
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3.  The  absolute  orientation. 


The  model  coordinates  can  be  transformed  to  the  coordinate 
system  of  the  object  -with  the  aid  of  one  scale  change,  three  trans¬ 
lations,  and  three  rotations.  For  this  purpose,  at  least  seven 
parameters  must  be  given,  for  instance,  two  completely  determined 
control  points  in  the  coordinate  system  of  the  object  and  one  addi¬ 
tional  coordinate  (in  terrain  mapping,  one  elevation  point). 

Certain  conditions  concerning  the  choice  and  location  of  the 
control  data  must  be  noted. 

It  should  be  noted  that  the  model  coordinate  errors  due  to 
various  sources,  for  instance,  the  elements  of  the  relative  orien¬ 
tation,  will  become  compensated  in  the  control  points  but  that  the 
elements  of  the  absolute  orientation  will  become  affected  with  errors 
accordingly.  Therefore,  the  errors  will  show  up  in  other  transformed 
model  points  together  with  the  original  model  errors. 

In  case  of  redundant  control  data,  a  certain  principle  must  be 
used  for  the  distribution  of  the  discrepancies.  The  method  of  least 
squares  is  also,  in  this  case,  of  importance  for  a  unique  solution 
of  the  problem  under  well-defined  conditions  for  the  error  propagation 
and  usually  leads  to  simple  calculations. 
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